Mon. Not. R. Astron. Soc. 000, [l]{T7] (2009) Printed 11 March 2010 (MN WT^ style file v2.2) 



The Possible Impact of GRB Detector Thresholds on 
Cosmological Standard Candles 

A. Shahmoradi^* and R. J. Nemiroff^* 

^Department of Physics, Michigan Technological University, Houghton, MI 49931 



ABSTRACT 

GRB satellites are relatively incfBcicnt detectors of dim hard bursts because they 
trigger on photon counts, which are number-biased against hard photons. Therefore, 
for example, given two bursts of identical peak luminosity near the detection threshold, 
a dim soft burst will be preferentially detected over a dim hard burst. This detector 
bias can create or skew an apparent correlation where increasingly hard GRBs appear 
increasingly bright. Although such correlations may be obfuscated by a middle step 
where GRBs need to be bright enough to have their actual redshifts determined, 
it is found that the bias is generally pervasive. This result is derived here through 
simulations convolving a wide variety of possible GRB brightnesses and spectra with 
the BATSE Large Area Detectors (LAD) detection thresholds. The presented analyses 
indicate that the rest- frame i^Fi, spectrum peak energy of long-duration GRBs, Ep^i„t, 
is not a good cosmological standard candle without significant corrections for selection 
effects. Therefore, the appearance of Ep ^nt in seeming correlations such as the Amati 
(Eiso — Ep^int), Ghirlanda {E^ — Ep^int), and Lisa — Ep^int relations is statistically real 
but strongly influenced by so far uncalibrated GRB detector thresholds. 
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1 INTRODUCTION 

The existence of correlations among the spectral parameters 
of Long-duration Gamma-Ray Bursts (LGRBs) has been 
touted as providing clues to the underling physics of GRB 
prompt emission and making LGRBs a useful tool for prob- 
ing cosmology in the distant universe. Recently reported at- 
tempts to use these correlations to construct a GRB Hubble 
diagram include those by Cardone et al. (2009), Amati et 
al. (2008a), Basilakos & Perivolaropoulos (2008), Cuesta et 
al. (2008), Liang et al. (2008), Schaefer (2007) - hereafter 
S07 - & Firmani et al. (2006). The investigation of possi- 
ble correlations among the parameters of LGRBs, however, 
dates back to the BATSE era when the cosmological origins 
of LGRBs was not yet established. Specifically in an early 
effort, Lloyd, Petrosian, & Mallozzi (2000) did an analysis 
to determine the degree of correlation between the observer- 
frame vFjj spectrum peak energy (Ep^obs) and the bolometric 
fluence of bright BATSE LGRBs (Sboi) and to investigate 
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to what extent the correlation is either intrinsic or cosmo- 
logical. After accounting for the data truncation due to the 
detection threshold, they concluded that there is probably 
a significant correlation between the rest-frame i/F,^ spec- 
trum peak energy (Ep^int) and isotropic-equivalent radiated 
energy {Etso)- 

While there is still no unique and robust interpreta- 
tion of these results (e.g. Levinson & Eichler 2005; Rees 
& Meszaros 2005; Eichler & Levinson 2004), the discovery 
of some outliers to these relations (e.g. Urata et al. 2009; 
Sugita et al. 2009; Bellm et al. 2008; McBreen et al. 2008; 
Campana et al. 2007; Rizzuto et al. 2007; Gehrels et al. 2006; 
Vaghuan et al. 2006; Sazonov et al. 2004; Soderberg et al. 
2004; Ghirlanda et al. 2004a - hereafter G04a) has raised 
two possibilities; that these correlations belong to only a 
sub-population of LGRBs, or that they are an artifact of 
the GRB detection process. The latter idea is bolstered by 
recent reports from several independent groups (Butler et 
al. 2009; Bagoly et al. 2009; Bagoly et al. 2008; Nava et 
al. 2008, hereafter N08; Butler et al. 2007, hereafter B07). 
Moreover, Nakar & Piran (2005a) (hereafter NP05a), found 



© 2009 RAS 



2 A. Shahmoradi and R. J. Nemiroff 



that a significant number of BATSE LGRBs are inconsis- 
tent with two of these relations known as the 'Amati re- 
lation' (Amati 2006; Amati 2002) and the 'Ghirlanda rela- 
tion' (Ghirlanda et al. 2007 - hereafter G07; G04a) which 
relate Ep^int of LGRBs to their isotropic-equivalent radiated 
energy {Eiso) and the coUimation-corrected energy {E~f) re- 
spectively. Following their analysis, Band & Preece (2005), 
analyzed a large sample of BATSE LGRBs and found that 
about 88% and L6% of their sample were inconsistent with 
the Amati and Ghirlanda relations respectively. Kaneko et 
al. (2006) - hereafter K06 - also reported an inconsistency 
of bright BATSE bursts with these correlations. 

Responding to these results, Ghirlanda et al. (2005a) 
argued that taking into account the intrinsic scatter of the 
Amati relation, the BATSE bursts may still be consistent. 
This claim, however, has also been challenged (Nakar & Pi- 
ran 2005b). Although all the above mentioned reports gen- 
erally conclude that the Amati and Ghirlanda relations are 
statistically non-compelling, the matter still remains con- 
troversial whether the reported correlations are completely 
due to selection effects in the detection process, or whether 
there is some real, statistically strong correlation between 
the LGRBs spectral parameters. It is noteworthy that Yo- 
netoku et al. (2004) - hereafter Y04 - have also tried to es- 
timate the redshifts of 745 BATSE LGRBs using a relation 
between Lisa and Ep^i„t- However, these estimates resulted 
in 21 GRBs being located beyond z > 12, and 35 even having 
no solution satisfying the Lisa — Ep^int relation, indicating 
that these 35 bursts show a large observer-frame peak en- 
ergy (Ep^obs) while being very dim. In addition, Tsutsui et 
al. (2008) have shown that the redshifts derived from the 
lag-luminosity & Lisa — Ep^i^t for 565 BATSE LGRBs are 
totally inconsistent with each other. 

Most recently, N08 and Ghirlanda et al. (2008) - here- 
after G08 - have reported the triggering threshold limits for 
several GRB detectors, including BATSE, SWIFT, Konus- 
Wind, BeppoSAX and HETE-II in the plane of peak energy 
vs. bolometric fluence [Sboi — Ep^obs) and peak energy vs. 
bolometric 1-second peak flux {Pboi — Ep^obs) of LGRBs. 
They also obtain the minimum fluence limits required for 
spectral analysis of a burst on these planes and conclude, 
contrary to the previous reports, that only 6% of BATSE 
LGRBs are certain outliers with respect to the Amati re- 
lation at > 3(T, while there are apparently no outliers to 
the Ghirlanda relation. They also find that the slope and 
the distribution of LGRBs are significantly different from 
the slopes of the curves of the minimum fiuences and peak 
fiuxes required for triggering and spectral analysis on these 
planes. Their analysis, however, is limited to about 380 
bright BATSE bursts analyzed by K06 and themselves which 
constitute only 14% of the entire number of GRBs detected 
by BATSE. The small sample of bursts used by N08 and G08 
makes it very hard to have an accurate investigation of the 
trigger threshold effects on the two planes of Sboi~Ep^obs and 
Pboi— Ep^obs- For this reason, they focus mainly on the limits 
imposed by the spectral analysis of LGRBs on these planes. 



leaving potential effects of trigger threshold limits on the 
distribution of bursts untreated. 

In this paper, while focusing our attention on BATSE 
LAD detectors, we study the trigger threshold effects on 
the joint distributions of the spectral parameters of BATSE 
GRBs using a method that avoids the difficult direct mea- 
surement of Ep^obs by substituting a much more easily found 
hardness ratio (Shahmoradi & Nemiroff 2010a). Based on 
Ep^obs estimates provided in Table (4) of Shahmoradi & Ne- 
miroff (2010a), we are then able to include 1900 BATSE 
GRBs, including ~ 500 short duration bursts (SGRBs) 
(Tbo < 3 [sec]) and more than 1400 LGRBs (Tgo > 3 [sec]) in 
our simulations of BATSE LAD detectors. In order to have 
an accurate investigation, we perform the analysis based on 
two entirely different approaches to GRB classification: 

(i) Based on the traditional GRB classification in which 
we label the bursts as SGRBs or LGRBs according to the 
duration division determined by Kouveliotou et al. (1993), 
but with a higher cutoff on the duration set at Tgo = 3 [sec] 
as compared to Tgo = 2 [sec] used therein. The higher cutoff 
is used to ensure minimal contamination of BATSE LGRBs 
by possible SGRBs with Tgo ~ 2 [sec] . 

(ii) Based on Fuzzy Clustering algorithms in which we as- 
sign each burst a probability of being LGRB or SGRB, de- 
rived from three distinct properties of Gamma-Ray Bursts: 
bolometric fluence (Sboi), spectral peak energy {Ep^obs) or 
equivalently, hardness ratio (HRh) as defined by Shah- 
moradi & Nemiroff (2010a), and the duration (Tgo). A de- 
tailed description of the methodology and algorithm used for 
the fuzzy classification of BATSE GRBs is given by Shah- 
moradi (2010). 

We are then able to proceed, in simulations and subse- 
quent analyses, to investigate the significance of the effects 
of triggering threshold limits on the detection and distribu- 
tion of BATSE GRBs in the planes of Sboi — Ep^obs as well 
as Pboi — Ep^obs and its possible effects on the Amati and 
Ghirlanda relations. Specifically, we show that the current 
realizations of the Amati relation as a low-dispersion {a ~ 
0.2 dex) Log-Log relation given by Amati (2006) & G07, are 
likely skewed by a large population of dim hard LGRBs that 
are missed due to the impossibility of a spectral analysis 
that determines Ep^obs- The significance of the detector 
threshold effects on these relations will be discussed in a 
separate paper by Shahmoradi & Nemiroff (2010b). 

Among all the gamma-ray observatories that have 
detected GRBs, BATSE provides the largest GRB database 
from a single experiment, consisting of observational data 
for 2704 GRBs. For many of the BATSE GRBs, high time 
and energy resolution data are available. The BATSE data 
are, therefore, the most suitable both in quantity and 
quality for detailed spectral studies and determination of 
the distribution of LGRBs on the planes of Sboi — Ep^obs and 
Pboi — Ep^obs, which can help to examine the veracity of the 
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reported correlations. 

The plan of the paper is as follows: i|2] is spent on the 
simulations of the BATSE Large Area Detectors (LAD) trig- 
gering thresholds and its possible effects on the bivariate 
distributions of GRBs in the planes of Stoi — Ep^obs and 
Pboi — Ep^obs- The possible reasons for the discrepancies be- 
tween our results and the findings of N08 and G08 will also 
be given. We discuss the tightness of the Ghirlanda relation 
and its connection with the Amati relation in !|3] The prop- 
erties of Amati relation in the observer and rest frames are 
investigated in i|4] The results will be summarized in i|5] 



2 BATSE DETECTION THRESHOLDS 
2.1 BATSE Trigger Algorithm 

The investigation of the BATSE triggering threshold, 
its specifications as well as its possible effects on the 
detection of different types of GRBs, is well documented 
in articles by the BATSE team (e.g. Pendleton et al. 1995; 
Pendleton, Hakkila & Meegan 1998; Paciesas et al. 1999; 
Meegan et al. 2000; Johnson, Meegan & Hakkila 2000; 
Brainerd et al. 2000; Hakkila et al. 2003a; Hakkila et al. 
2003b). Model-dependent studies of the BATSE triggering 
threshold for different Ep^obs - based on the observational 
data as well as comparison with other detectors - has 
also been presented by Band (2006), Band (2004), Band 
(2003) & Band (2002). A precise determination of the 
trigger threshold, however, involves modeling the triggering 
algorithm of BATSE and its observational efficiency, taking 
into account the sensitivity of BATSE LADs at different 
angles of the incident photons. 

BATSE had a relatively simple triggering algorithm 
compared to other instruments such as BAT onboard 
SWIFT and HETE-II (e.g. Band 2003). It generally 
triggered on a burst if the count rate in the second most 
brightly illuminated detector exceeded a threshold specified 
in units of standard deviations - nominally 5.5a for normal 
incidence on a single detector - above background rate 
that was determined for each detector over a given time 
interval, usually set at 17.4 sec. However, as indicated by 
Band (2003), the requirement that a trigger occurs when 
the fiux in the second most brightly illuminated detector 
reaches 5.5a, raises the significance threshold from 5.5a 
- for a normal incidence - to 5.96a — 7.78a above the 
background, depending on the cosine of the angle between 
the normal axis of the detector and the source. BATSE 
was programmed to trigger on any of three time scales: 64 
ms, 256 ms and 1024 ms and the trigger energy range was 
generally set to 50 KeV - 300 KeV. 

To simulate BATSE, we use BATSE's Detector Re- 
sponse Matrices (DRMs) and public GRB data available 



through the online HEASARC archives The fluences for 
four BATSE channels covering the energy range of 20 — 2000 
KeV - which we treat as bolometric - were taken from the 
"Current BATSE" catalog. Jimenez et al. (2001) have shown 
that the fluences resulting from the processing pipeline used 
to create the BATSE catalog are slightly - up to a factor of 2 
- larger than the fluences obtained by fltting high-resolution 
spectra. By simulation, using a wide range of Ep,obs (1 — 2000 
KeV) with typical photon indices a = —1.1 and /3 = —2.3, 
we find the ratio of the bolometric fluence (given the rest- 
frame energy range of 1 — 10, 000 KeV) as defined in G07 
and the fluence in BATSE energy range to be very close to 
unity (on average < 1.2). Comparing this ratio with the re- 
sults of Jimenez et al. (2001), extension of the energy band 
to compute bolometric fluence via spectral fits appears to be 
unnecessary for the current analysis (Band & Preece 2005; 
Friedman & Bloom 2005). 

In some previous analyses (e.g. N08 & G08), the trig- 
ger threshold limits were generally obtained using the Band 
model with typical photon indices fixed to a — —1.1 & 
P = —2.3 as a representation of the spectra of the whole 
sample of BATSE LGRBs. However, we will show in §2.2| 
& j ]2.3| that these types of simulations are strongly affected 
by model biases and circularity problems. Instead, we rely 
on the light curves of the GRBs in their original forms as 
detected by BATSE detectors and use them to find the 
lower limits for the trigger threshold on Sboi — Ep^obs & 
Pboi —Ep^obs planes by decreasing monotonically the amount 
of received photon counts in all energy channels for each in- 
dividual burst. 

The analysis of background count rates in the light 
curves of 2145 BATSE GRBs gives average minimum 
required peak fluxes of 0.66-0.86, 0.34-0.43, 0.18-0.24 
ph cm~^ for the three timescales 64, 256 and 1024 ms 
respectively, given the aforementioned BATSE significance 
threshold. These values correspond well to the lowest ob- 
served peak fluxes of BATSE GRBs on these timescales. We 
also include all three timescales in the trigger criteria of our 
simulation. This is particularly important for SGRBs, as 
they are mainly triggered on 64 ms and 256 ms timescales. 
The entire simulation algorithms in this work are written in 
FORTRAN. 



2.2 Sboi - Ep.obs Plane of BATSE Bursts 

Figure [l] shows the Sboi— Ep^obs plane of 1900 BATSE bursts 
with continuous light curves and known durations and flu- 
ences, chosen from the complete BATSE catalog. Interest- 
ingly, the distributions of both SGRBs (blue dots) & LGRBs 
(red dots) look similar to each other and exhibit similar 
slopes on both their left and right sides (the dotted lines 
in Figure [T}. The only difference is that SGRBs are shifted 
towards the upper left of the Sboi — Ep^obs plane, possibly 

^ http:/ /gammaray.msfc. nasa.gov/batse/grb/catalog/current/ 
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Figure 1. Plot of the spectral peak energy, Ep abat versus the 
bolometric fluence, Sf,oi, for 1900 BATSE GRBs with fluences, 
duration, and light curves available in the BATSE GRB cata- 
log. The red and blue dots represent LGRBs (Tgo > 3 [sec]) and 
SGRBs (Tgo < 3 [sec]) respectively. The dotted lines are the ap- 
proximate left and right boundaries for the distribution of each 
population (SGRBs & LGRBs) obtained via a method described 
in §2.2[ Both SGRBs and LGRBs appear similarly bounded on 
both the upper left and lower right. 



indicating that the duration of the burst has an important 
effect on its detection and determines the position of the 
GRB on the plots of Ep^obs versus the bolometric fluence 
(Sboi)- Moreover, the scatter of SGRBs on this plot appears 
to be smaller than the scatter of LGRBs. SGRBs also show 
a slightly tighter correlation (Kendall tk ~ 0.43, ISa) com- 
pared to LGRBs (tk — 0.32, 18a). This could be due to the 
narrower duration distribution of SGRBs as compared to 
LGRBs. 

Since the duration of the burst plays an important role 
in its detection and also in connecting the triggering peak 
flux of the burst to its fluence, the effects of GRB duration 
should also be considered in the studies of selection effects 
on Sboi — Ep^obs plane of GRBs. To expand on this, we have 
created a simulation where the BATSE GRB sample is di- 
vided into eight subgroups with respect to their Tgo dura- 
tions. This simulation takes a given GRB and, while keeping 
duration constant, artificially decreases its fiuence uniformly 
across the light curve until the GRB would no longer trigger 
BATSE. The results of this simulation are shown in Figure[2] 
Each gray star shows the lowest fluence point for one partic- 
ular BATSE GRB in the sample. Any point to the right of a 
gray star is detectable by BATSE if the burst has the same 
spectral parameters of the gray star, but a higher fluence. 

One noticeable feature of the first four graphs on 
LGRBs in Figure [2] is that the slopes of the boundaries on 
the left and right sides of each GRB subsample appear dif- 
ferent. The sense of this difference is that the slope on the 



left side appears significantly steeper than the slope on the 
right side. The slopes at the borders of the distributions are 
obtained via an iterative algorithm that looks for the best 
fit line to the border which minimizes the correlation of the 
residuals of the boundary points to that line with the dis- 
tances of their projections on that line from a conventional 
point on the same line chosen as the origin for calculating 
the distances. The algorithm is repeated until the correla- 
tion coefficient reaches a level that the two parameters can 
be regarded as totally independent variables. 

The difference between the left and right slopes for the 
LGRB plots in Figure [2] has an important implication. In 
particular the left slope might be attributable to trigger 
threshold limits. It is generally argued (e.g. N08 & G08) 
that the existence of a strong correlation in the plane of 
Sboi — Epochs is an indication of even stronger correlation 



in the plane of Ei, 



Although this is true, the ex- 



istence of an intrinsically tight. Log- Log linear correlation 
between Ep^tnt & Eiso must essentially result in an observed 
distribution of LGRBs on Sboi — Ep^ots plane that has nearly 
the same slopes on its left & right boundaries (e.g. Figure 5 
& 6 of Bosnjak et al. 2008). This is contrary to what is seen 
in the plots of Figure [l] and Figure [2] 

Another important attribute of the plots in Figure [2] is 
the negative correlation between the duration of the bursts 
and the fraction of 'certain' outliers to the Amati relation at 
> 3cr. The outliers are determined following the method de- 
scribed by NP05a and Band & Preece (2005), also used ear- 
lier than these authors by Ghirlanda et al. 2004b for SGRBs. 
This negative correlation has possibly no physical origin and 
is due to selection effects caused by the durations of the 
bursts, that is, longer duration bursts are generally less de- 
tectable than shorter duration bursts of the same fluence. 
In other words, given the same fluence for two GRBs, the 
longer duration burst will have a larger Fluence to Peak flux 
Ratio (FPR) than the shorter duration burst which makes 
it less detectable for GRB detectors. 

The results of our simulations are clearly different from 
the findings of N08 & G08 who obtain a distribution and 
slope for their sample of LGRBs that are far from BATSE 
trigger limits. The possible reasons for such discrepancy are 
summarized below: 

(i) Circular Logic Problem: In the method used by 
N08 & G08, the authors rely on the data from a fraction 
of LGRBs ^detected' by BATSE and 'spectrally analyzed' 
to constrain the parameters involved in their simulations. 
These parameters, such as FPR, are required to relate the 
peak flux of the simulated burst to its fiuence, and the du- 
ration of the burst. 

The use of the already detected GRBs to obtain the lim- 
iting parameters for their simulations, however, causes their 
analysis to suffer from a circular logic problem. Specifically, 
the value of FPR that they use (~ 6), is representative of 
detected GRBs, not the entire population of GRBs which in- 
cludes both detected and undetected bursts. Figure [4] shows 
a frequency histogram of FPR values for 1053 bright BATSE 
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Figure 2. Plots of Ep ^bs versus S'^oj for different duration subsets of real BATSE (black points) and simulated BATSE (gray star) data. 
Each simulated gray star represents the theoretically dimmest visible point for a particular BATSE GRB, below which BATSE would 
not have triggered on the burst, given BATSE detector thresholds. The dashed-red line in each plot represents the 3o" upper limit to the 
Amati relation (see §3.2), above which any point is a certain > 'ia outlier to the Amati relation given by G07. 
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Tgo [sec] 

Figure 3. Plot of the fraction of > Scr outliers to the Amati 
relation versus Tgo duration for 1900 BATSE GRBs, represented 
by the hatched histogram. A negative correlation is apparent, 
indicating that regardless of the type of the burst, whether LGRB 
or SGRB, the number of outliers to the Amati relation increases as 
the durations of the bursts decrease. The solid line represents the 
normalized histogram of Tgo of the LGRBs used by G08 to define 
the low-dispersion Log-Log linear Amati relation. As clear in the 
graph, most of the LGRBs in the G08 sample are chosen from a 
specific subsample of LGRBs that have very long durations. The 
background (gray) histogram represents the relative frequency 
histogram of Tgo for 2041 BATSE GRBs. Without recourse to any 
statistical significance tests, it is clear that the sample of GRBs 
used to define the Amati relation does not represent the 
entire population of Long-duration GRBs and it is biased 
towards the longest duration GRBs. 




FPR 

Figure 4. Histogram of FPRs, Fluence to 1-sec Peak fiux 
Ratios, for BATSE LGRBs. Bolometric FPRs are plotted as the 
gray histogram, while the solid-line histogram includes FPRs cal- 
culated in the BATSE detection energy range 50-300 KeV. Both 
histograms peak at FPR 16, in contrast to the peak value 
reported by G08 {FPR 6) obtained for a smaller sample of 
BATSE LGRBs, represented by the vertical solid line in the 
graph. Only 27% of the whole sample have FPR < 6. Under- 
estimating this parameter can result in underestimation of the 
significance of BATSE detector thresholds on the distribution of 
LGRBs in Si,oi — £^p,o6s plane, as explained in |2.2| This likely 
explains one of the reasons for the significance of the differences 
between the results of the simulations of BATSE detector thresh- 
olds presented here (Figure [2]| and the results obtained by other 
authors (e.g. Figures 1, 2, 5 & 8 in N08 & Figure 5 in G08). 



LGRBs (Tgo > 3 [sec]) with nonzero fluences in all 4 en- 
ergy channels. Here FPR is defined as the ratio of fluence to 
1-second peak flux (Pio24). Both fluence and the peak flux 
were calculated first in 50 — 300 KeV energy range (the solid- 
line histogram of Figure |4| and next for the whole BATSE 
energy range (the gray background histogram). The graph 
indicates a peak at FPR ~ 16, a factor of 2.5 higher than 
the FPR assumed in the simulations of N08 & G08. The 
visible asymmetry in the histogram could be partly due to 
selection effects on high FPR LGRBs as they are generally 
less detectable compared to low FPR LGRBs. 

(ii) Duration Effects on GRB Detection: Another 
factor, possibly responsible for the discrepancies between 
our results and the results obtained by N08 & G09, is the 
range of durations of the bursts considered in the simula- 
tions. N08 & G08 use 5 seconds & 20 seconds as the lower 
and upper limits. Such an assumption, however, underesti- 
mates the effect of threshold limits on very long duration 
GRBs that constitute a non-negligible fraction of the entire 
sample of LGRBs. In particular, for the BATSE sample of 
LGRBs considered in this work, the total number of bursts 
with 3 [sec] < Tgo < 20 [sec] is only 33% of the whole. This 



could be even lower, considering the fact that very long dura- 
tion LGRBs have less chance to be detected. For LGRBs de- 
tected by the Burst Alert Telescope (BAT) onboard SWIFT 
satellite which is more sensitive to longest duration GRBs 
compared to BATSE LADs, this ratio decreases to 28% (Fig- 
ure |5|. 

In order to show the effects of duration on GRB trig- 
gering, we made a contour plot of Ep^obs — ^'1024 — Tko for 
all BATSE GRBs by normalizing their fluences (Sboi) to a 
canonical fluence level (1Q~^ erg cm~^s~^) in Figurelsl Here 
we use T50, since it is a more accurate measure of GRB du- 
ration than Tgo which is more dependent on the background 
fitting model and the choice of burst start and stop times, 
that are typically set by hand (e.g. Paciesas et al. 1999; 
B07). The graph clearly shows that for the same fluence 
and Ep^obs, longer duration bursts have lower peak fluxes 
than the shorter bursts. In other words, given the same flu- 
ence & duration for two bursts, the harder burst will have 
a lower chance of a triggered detection than the softer one. 
Were the above points considered in the analyses of N08 
and G08, their Trigger and Spectral Threshold (TT & ST) 
curves for BATSE would shift towards the distribution of 
BATSE LGRBs in then plots of Sboi - Ep,obs (e.g. Figure|8|. 



© 2009 RAS, MNRAS 000,[T]{T7] 



The Possible Impact of GRB Detector Thresholds on Cosmological Standard Candles 7 



0.10 



■D 
03 
N 

re 
E 



3s < Tgp < 20s 

BATSE: 33% 
SWIFT: 28% 
G08: 25% 





'90 



[sec] 



-p.obs 



[KeV] 



Figure 5. A histogram of Tgo durations for 1390 BATSE LGRBs 
with Tgo > 3s (gray background histogram) as well as 345 SWIFT 
LGRBs with Tgo > 3s (hatched foreground histogram). Also 
shown is the Tgo histogram of 76 LGRBs used in G08 to con- 
struct the latest update of Amati relation, here normalized to 0.1 
(rather than 1) for a better comparison of the three histograms. 
The solid vertical line represents Tgo = 20s which is used by G08 
&; N08 as the upper limit for the durations of the bursts in their 
simulations of Trigger and Spectral analysis Threshold (TT & 
ST) hmits. However, only 33% of the entire sample of BATSE 
LGRBs have Tgo < 20s. This ratio reduces to 28% for the cur- 
rent SWIFT sample of LGRBs and to 25% for G08 76 LGRBs 
used to construct the Amati relation. The latter is particularly 
important in the analyses done by G08, since they compare TT 
& ST limits of different instruments obtained for Tgo < 20s with 
a sample of LGRBs that 75% of them have Tgo > 20s that results 
in a significant underestimation of the selection effects. 



The two parameters - FPR & duration - however, do not 
significantly affect the shapes of TT & ST limits, since they 
only result in a shift in the positions of the entire TT & ST 
limits on the plane of Sboi — Ep^obs- The important factors in 
determining the shapes of these limits are discussed in §2.3| 



It is notable that two famous outliers of the Amati re- 
lation, GRB 980425 and GRB 031203, both lie in a region 
where the sample of G08 LGRBs exist and therefore cannot 
be flagged as certain > 3(j outliers to the Amati relation 
(G07) via the method given by NP05a and Band & Preece 
(2005). This is a major weakness of the method given by 
these authors as indicated by themselves, since it can only 
set a lower limit for the number of certain outliers to the 
Amati and Ghirlanda relations. Nevertheless, following this 
method, for any linear relation among the rest-frame spec- 
tral parameters of GRBs - Xrest & Yrest - we can write. 



Log{Yrest) = Q + l3Log{Xrest) + ^cr. 



(1) 



where ^ stands for the significance level of being an outlier 
to a particular relation considered in the method. Here, we 
use ^ = 3, corresponding to Scr level (e.g. Figures [2]&[7]l, also 



Figure 6. A plot of the observed peak flux, Pio24, versus the 
spectral peak energy, Ep obat for 1900 BATSE GRBs segregated 
in T50 by color. Each color region represents a population of the 
bursts with similar T50 durations. The fluences of all bursts are 
normalized to a fiducial fluence level (10~^ erg cm~^s~^). The 
vertical axis is the 1-scc peak flux of GRBs (P1024) and the dot- 
ted line represents the nominal BATSE trigger threshold on 1024 
[sec] timescale. For the same fluence and hardness (either Ep ^bs 
or HRh), bursts with longer durations have peak fluxes that arc 
orders of magnitude less than the peak fluxes of shorter duration 
bursts, and therefore will be less detectable. This further illus- 
trates the selection effect bias discussed in 32.21 



a range of values for ^ as depicted in Figure [7] According 
to the assumed significance level 5, Eqn. ^ can then be 
transformed into the observer frame, knowing that. 



Yrest 
Yobs 

X 
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where friz) & .fx{z) are functions of redshift that relate 
the observer-frame spectral parameter to the corresponding 
rest-frame parameter. For the Amati, Ghirlanda, and Yone- 
toku relations, Eqn. ||3| can be written as, 
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where is the beaming factor which is determined ob- 
servationally from modehng the evolution of the afterglow 
for each individual GRB, di, is the luminosity distance in 
concordance cosmology assuming = 0.27, Ha ~ 0.73, 
Ho = 72 Km/sMpc (Komatsu et al. 2010) & c as the speed 
of light. 

For the Amati and Ghirlanda relations, the function 
A{z, a, P, 5, a) has a maximum, Amaxict, /?, ^, a). Any burst 
with spectral parameters such that, 



xt 



(9) 



will be a certain outlier to these relations at > ^ct for any 
redshift it might have. Using the above method we find that 
at least 19% (or 21%) of BATSE LGRBs are certain 
outliers to the Amati relation as given by Ghirlanda 
et al. (2008) - (G08) - at > 3a significance level, 
based on the classical definition — Tgo > 3 [sec] — (or 
fuzzy clustering classification) of BATSE LGRBs. 
Moreover, the consistency of the BATSE LGRBs 
with the Amati relation of G08 is strongly rejected 
with KS significance probability of 10"^^" (Figure [t]). 
This fraction could be even higher knowing that the current 
sample of 1900 GRBs used here is only 2/3 of the total num- 
ber of GRBs detected by BATSE, while the remaining GRBs 
- not presented here - are generally bursts with very low flu- 
ence close to trigger threshold limits, excluding exceptional 
bursts such as those with data gaps in their light curves. 

In addition, it can be shown that the Amati relation 
as given by G08 is heavy on the side corresponding to the 
region of dim, hard bursts. Assuming a Gaussian distribu- 
tion of the data around the best linear fit as considered by 
G07, one would expect to observe 50%, 16% and 0.023% of 
the bursts lying at > ^cr for ^ = 0, 1, 2 respectively, on the 
dim-hard side of the rest-frame Amati relation. However, 
the observed fractions are 80%, 60% & 37%, resulting in an 
observed excess of the bursts > 30%, > 44% & > 35% at 
> ^(T, ^ = 0, 1,2 respectively on the dim-hard side of the 
Amati relation. The expected & observed fractions are plot- 
ted for a continuous range of significance levels in Figure [7] 
{Lower Left & Lower Right). 



2.3 Pboi - Ep^obs Plane of BATSE Bursts 

Several authors have reported a strong correlation between 
Ep^int and the rest-frame isotropic peak luminosity (Lisa) 
of LGRBs (e.g. Ghirlanda et al. (2009), hereafter G09; S07; 
Ghirlanda et al. 2005b; Y04; Schaefer 2003). The relation 
was originally constructed from a handful of BeppoSAX and 
BATSE bursts with known redshift and was updated by S07 
for a sample of 64 LGRBs detected mainly by SWIFT and 
HETE-II and most recently by G09. The study of selection 
effects for this relation however, is much older than the re- 
lation itself (e.g. Lee & Petrosian 1996; Lloyd & Petrosian, 
1999; Lloyd, Petrosian & Mallozzi 2000). 

Following the same method used for investigating the 



selection effects in the plane of Stoi — Ep^obs in |2.2[ we 
use the light curves of BATSE LGRBs and their observed 
1-second peak fluxes to study the bivariate distribution of 
GRBs in the Pboi — Ep^oba plane. This is done by gradually 
decreasing the normalization constant of the spectrum for 
each LGRB up to the limit where BATSE could not trigger 
the burst, based on the minimum required peak fluxes on the 
three different timescales that BATSE used for triggering: 
64, 256 and 1024 ms. The results are shown in Figure[7]( Top 
Right) . 

Similar to the Sboi — Ep^obs plot (Figure [2|, inspection 
of Figure [7] indicates that the far left side of the distribution 
of BATSE LGRBs on the plot of Pboi - Ep^obs is affected 
by the trigger threshold. Additionally, there appears to be 
a lower limit for Ep^oba, given an observed fluence, on the 
far right side of the sample. Since maximum brightness is 
not limited at detection, this is unlikely to be due to data 



truncation (also depicted in Figure 10 1. Moreover, the left 
and right boundaries of the sample have different slopes: 
mi = 1.57 ± 0.05 for Ep^obs < 400 KeV & m, = 0.89 ± 0.14 
for Ep^obs > 400 KeV, = 0.62 ± 0.02. 

The GRB samples used by S07 & G09 to construct the 
Liso ~ Ep^int relation, have correlation coefficients of tk ~ 
0.48 ± 0.07(5.4cr) k tk = 0.48 ± 0.05(6. 7a) respectively in 
the observer frame, while the sample of 1053 BATSE LGRBs 
considered here has tk = 0.36 ± 0.02, 17cr. Following the 
same procedure as used for the Amati relation, we find that 
at least 8% of BATSE LGRBs are outliers to Li^o - 
Ep,int of GOO at > 3(T level for both the traditional 
{Tgo > 3 [sec]) or the fuzzy classification of BATSE 
LGRBs. In addition, the consistency of the BATSE 
LGRBs with Liso — Ep.int of G09 is strongly rejected 
with KS significance probability of 10^^*" (Figure [t]). 



It can be shown that the Li, 



-^Pyint 



relation is heavy 



on the side corresponding to the region of dim, hard bursts. 
Assuming a Gaussian distribution of the data around the 
best linear fit as considered by N08, one would expect to 
observe 50%, 16% and 0.023% of the bursts lying at > ^u, 
^ — 0,1,2 respectively, on the left side of the rest-frame 
Liao — Ep^int relation. However, the observed fractions in- 
dicate an excess of the bursts > 40%, > 43% & > 21% 
at > ^(T, ^ = 0,1,2 respectively on the dim-hard side the 
Liao ~ Ep,int relation. 

The results of the simulations shown in Figure [t] ( Top 
Right) again indicate the contamination of BATSE LGRBs 
distribution on Pboi—Ep,obs plane by trigger threshold limits. 
This is again contrary to the findings of previous authors, 
in particular N08 who find that the trigger threshold (TT) 
and spectral threshold (ST) limit curves for BATSE, as well 
as other instruments, are quite far from the bivariate distri- 
bution of the LGRBs on Sboi — Ep^oba & Pboi — Ep^oba plots. 
The reasons for the discrepancies are likely as follows: 

In ^ model- dependent' studies of selection effects it is usu- 
ally assumed that all types of LGRBs with different Ep^oba, 
have the same low and high energy photon indices on av- 
erage, typically a ~ —1.1 & /3 ~ —2.3. In these studies, 
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Figure 7. Top Left: Data plot of Ep^obs versus Stoi for 1426 BATSE LGRBs among all 1900 BATSE GRBs, classified via a fuzzy 
clustering algorithm based on their bolometric fluence, hardness and duration. BATSE LGRBs are represented by the green dots. Also 
shown are the LGRBs of G08 used therein to construct the Amati relation, represented by the red & blue dots, corresponding to LGRBs 
detected by SWIFT and other instruments, respectively. The solid line represents the 3(t upper limit for the Amati relation as reported 
by G07. The 3 black stars represent GRB 980425 & GRB 031203 (outliers to both the Amati and Ghirlanda relations), also GRB 070125, 
consistent with the Amati relation but an outlier to the Ghirlanda relation at 5(t, as reported by Bellm et al. 2008. The black filled 
triangle represents GRB 060505, a far outlier to Amati relation. The gray star represents the sub-luminous SWIFT GRB 050826 (B07), 
the only LGRB - with measured redshift - that can be flagged as an outlier to the Amati relation via the method given by NP05a, 
regardless of its redshift. The white head-down triangle represents XRF 050416A that Sakamoto et al. (2006) reports its consistency with 
the Amati relation as a further confirmation of the physical origin of this relation. We show in j|4]that as long as any burst - including 
XRF 050416A - is chosen in the region where the current sample of LGRBs defining the Amati relation reside in the observer frame, 
it will be consistent with the Amati relation at < 3ct for any redshift z > 0.09. The red star represents GRB 060218 for the spectral 
parameters given by B07. Top Right: A plot of Ep^obs versus Pf,oi for 1053 bright BATSE LGRBs with cataloged nonzero fiuence in 
all BATSE energy channels. The black dots represent real BATSE GRB data, while the red dots represent the same BATSE bursts 



artificially reduced in fluence to the BATSE trigger threshold level. The solid line is the 3cr upper limit for Li, 



relation given 



by Ghirlanda et al. (2009). While there is seemingly no selection effect present on the right side, the far left side of the sample appears 
to be affected by the BATSE trigger threshold limits. Bottom Left: Plot of the fraction of BATSE LGRBs that are certain outliers 



to the Amati & Li. 



relations — represented by the red & blue curves respectively — at the given significance levels. The solid 



curves are obtained based on the classical definition of LGRBs (Tgo > 3 [sec]) while the dotted curves are based on the fuzzy clustering 
classification of 1900 BATSE GRBs (Shahmoradi 2010). The green curve represents the expected fractions for both relations assuming 



Gaussian residuals for both the Amati relation given by G08 & Li, 



relation given by G09. Bottom Right: The absolute 



differences between the observed and the expected fractions of BATSE LGRBs (Tgo > 3 [sec]) - represented by the solid rod, blue & 
green curves in the ^ Lower Left' plot - at the given significance levels. Based on the traditional (or fuzzy clustering) classification 
of BATSE GRBs, at least 19% (or 21%) & 8% (or 8%) of the BATSE sample of LGRBs considered in this work can 



be flagged as certain outli 
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^p,int relations respectively at > 3(T significance level. According 
) test, the similarity of the expected and the observed fractions for both the Amati & 
relations is rejected at high significance levels corresponding to KS significance probabilities of 10"^^" &c 
respectively. Such extremely small probabilities invalidate the common assumption of Gaussian distribution of the residuals 
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scatters for these relation than those proposed by Amati et al. (2006), G08 &i G09. 
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the trigger as well as the spectral analysis threshold limits 
in Sboi — Ep^obs & Pboi — Ep,obs planes are determined via 
a presumption that all LGRBs can be well described by a 
typical spectral model - usually the Band model - such as 
those done by G08 & N08. This is true when the high- and 
low-energy photon indices are statistically independent of 
Ep^obs- G08 find no dependency of a to Ep^obs for a sample 
of LGRBs used to construct the Amati relation therein. 

However, Shahmoradi & Nemiroff (20i0a) have shown 
that there are likely significant - and in some cases strong - 
positive correlations among the high- & low-energy photon 
indices of the three GRB models: Band, COMP (GPL) & 
SBPL and Ep^obs of bright BATSE GRBs. 

These positive correlations among the spectral param- 
eters can have significant effects in model- dependent studies 
of selection effects. To illustrate this, we also simulated the 
BATSE trigger limits following the model-dependent meth- 
ods as presented by G09, G08, N08 & Band (2003). Figurejs] 
[Right) is a plot of Trigger Threshold (TT) curves assum- 
ing three different models (Band, SBPL & GOMP) for the 
GRBs with fixed photon indices, given in the caption of the 
figure. 

The use of fixed averaged photon indices for all bursts 
with different Ep^obs clearly results in an underestimation of 
the minimum required bolometric peak flux {Pboi) to trigger 
a bursts at high Ep^obs- If this point is considered in ^model- 
dependent^ simulations of BATSE TT and ST curves, the 
slopes of GOB & N08 limiting curves would show much more 
proximity to the distributions of the bursts on Pboi — Ep^obs & 
Sboi — Ep^obs plots. Among all the TT curves obtained for 
different spectral models, the GOMP model curve with a — 
—0.5 (black dashed line in the Right plot of Figure |8| best 
outlines the dimmest points of Pboi — Ep_obs plane detectable 
by BATSE (green stars in the same figure). 



3 THE REALITY OF THE GHIRLANDA 
RELATION 

Although the present and the previous analyses done by sev- 
eral authors (e.g. Butler et al. 2009; B07; NP05a; Band & 
Preece 2005) strongly suggest that the Amati relation is due 
to complex selection effects not only in triggering process 
but also in the spectral analysis and redshift measurement, 
these cannot rule out other explanations for the inconsis- 
tencies. Several alternative explanations were discussed by 
Ghisellini et al. (2006), the most popular being the off-axis 
model for GRBs that appear to be outliers to Amati rela- 
tion. Therefore, as Amati (2008b) has suggested, the exis- 
tence of low-dispersion - as compared to Amati relation - 
three-parameter correlations, such as the Ghirlanda relation 
and the empirical relation Ep^i„t — Eiso — tjet given by Liang 
& Zhang (2005), could strongly favor a physical origin of the 
Amati relation and possibly other GRB correlations. To in- 
vestigate such a possibility, we first determine the number 
of certain outliers to the Ghirlanda relation for the homo- 
geneous medium case as given by G07 under the extreme 



condition that all GRBs have a beaming factor fb — 1. Fol- 
lowing the same method used for the Amati & Lisa — Ep^int 
relations in the previous sections, we find only <0.6% of the 
whole BATSE sample of LGRBs presented here to be certain 
outliers to the Ghirlanda relation. The fraction of outliers is 
therefore a significant difference between the Ghirlanda and 
Amati relations. In this section, we look into this difference 
more deeply. 

The Ghirlanda relation was initially presented by G04a 
for a sample of 16 LGRBs with jet opening angles esti- 
mated from the achromatic break of their afterglow light 
curves. They reported a significant scatter reduction and 
correlation improvement when the data were transformed 
form Eiso — Ep^int to E-y — Ep^int plane, with Spearman's 
rank correlation coefficients of = 0.80, — 0.94 and 
dispersions of 0.57 dex and < 0.1 dex about the best fit 
lines of the Amati and Ghirlanda relations respectively. In- 
terestingly, Ep^int is highly correlated with both Eiso & E^ 
in G04a sample. However, the hope for the existence of a 
significant difference between the correlation coefficients of 
the two relations vanishes when it is found that the Am- 
ati and Ghirlanda relations, considering the same sample 
of bursts (i.e. 16 GRBs with measured jet opening angle 
given by G04cl) for both relations, have correlation coeffi- 
cients that are within the la uncertainties of each other. 
The uncertainties in the correlation coefficients were de- 
termined via the bootstrap method by generating a large 
enough number of synthetic data sets. The resulting Kendall 
rank correlation coefficients were tk,a = 0.65 ± 0.16 (3.5(j) 
& TK,G = 0.85 ± 0.08 (4.6cr) for the Amati & Ghirlanda 
relations respectively. Moreover, the scatter of the Amati 
relation reduces to a a ~ 0.14 when fit to the same sample 
of 16 LGRBs, which is comparable to the ac ~ 0.08 for 
the Ghirlanda relation . The apparent correlation difference 
between the two relations diminishes yet further when con- 
sidering the whole sample of GRBs given in G04a, for which 
TK,A = 0.74±0.08 (5.2cr) & tk.g = 0.80±0.07 (5.6o-). Both 
relations have the same scatter (cr = 0.15) about their best 
linear fits. GRB 970508 was excluded from the above analy- 
sis because of its uncertain Ep^obs in G04a ranging from 145 
KeV to >800 KeV, 

The small size of the G04a sample with firmly reported 
9j^t itself raises questions about the correlation improve- 
ment of the Ghirlanda relation. Furthermore, were the cited 
correlation improvements to have a physical origin, they 
should manifest themselves more strongly in larger samples 
of GRBs. Unfortunately, a significantly larger sample is not 
yet available. 

A recent update of the Amati, Ghirlanda & Liang- 
Zhang relations have been given by G07, extending the num- 
ber of GRBs with firmly measured spectral data from 16 
in G04a to 24 in G07. Reanalyzing the sample of LGRBs 
given in G07, we confirm the correlations and scatters found 
therein: cja = 0.20, aa = 0.09, glz = 0.10 for the Amati, 
Ghirlanda & Liang-Zhang respectively. However, consider- 
ing the same sample of bursts (i.e. only those with firmly 
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Figure 8. Comparison of the simulated BATSE trigger limits based on the BATSE detected GRBs and the simulated 
'spectral-model dependent' BATSE trigger limits of G08 and N08. Left: Plot of Ep ^bs versus Stoi demonstrating trigger 
threshold effects for BATSE. The colorful stars represent actual BATSE LGRBs (T90 > 3 [sec]) that are artificially dimmed in fluence 
to the detector threshold and so best show the real sensitivity of BATSE to dim LGRBs. By contrast, the two superposed curves - 
connected to each other by horizontal lines - show the effective Trigger Threshold (TT) limits obtained following G08 who use the Band 
model with fixed average photon indices, a = —1.1 & = —2.3, and FPR = 6 in their analyses to estimate BATSE triggering threshold. 
The two curves correspond to the TT limits for Tgo = 5 [sec] (left) & T90 = 20 [sec] (right) as given in Figure 5 of G08. The plot indicates 
that G08 significantly underestimates the effects of BATSE trigger threshold on the distribution of BATSE LGRBs in Si,oi — E^y obs plane 
(ij2.2|l. Right: Plot of Ep ^bs versus P(,oj demonstrating trigger threshold effects for BATSE. The green stars represent 1053 bright 
BATSE LGRBs that are artificially dimmed in fluence to the detector threshold and so best show the real sensitivity of BATSE to dim 
LGRBs. By contrast, the superposed lines show the effective trigger threshold limits obtained using different spectral models. The solid 
lines represent the TT limits for average photon indices of K06 sample of LGRBs with a = —1.1 h fi = —2.3, Ai = —1.3 & A2 = —2.3, 
a = —1.3 for three spectral models Band, SBPL &: COMP respectively. The dashed lines are TT limits obtained using higher photon 
indices with a = —0.5 & /3 = —2.0, Ai = —0.5 & A2 = —2.0, a = —0.5 in the same order as mentioned. For the SBPL model, the fifth 
parameter (break scale) was fixed to the average value of K06 GRB sample (A = 0.2) in both cases. None of the effective TT lines closely 
follow the actual TT points. 



measured spectral data, including 9j^t) for all three rela- 
tions, the scatter of Amati relation becomes comparable to 
the two others [a a = 0.14). Also, the correlation coefficient 
improvement observed in 16 GRBs of G04a sample deteri- 
orates substantially: tk.a = 0.76 ± 0.09 (5.2a") & tk,g ~ 
0.82 ± 0.06 (5.6cr). Including GRB 070125, GRB 071010B 
& GRB 050904, recently found outliers to the Ghirlanda re- 
lation at > So- level (Urata et al. 2009; Sugita et al. 2009; 
Bellm et al. 2008), makes the Ghirlanda relation comparable 
to the Amati relation and results in tk,a = 0.76±0.08 (5.5cr) 
& TK,G = 0.73 ± 0.07 (5.4a) with aA = 0.16, aa = 0.22. 

It is important to mention that we did not include any 
of the G07 bursts with uncertain spectral parameters (such 
as 9 jet or Ep^obs) in the above analysis. The inclusion of 
all 33 bursts (by fixing the parameters with lower or upper 
limits to the values given) in G07, however, results in an 
even higher dispersion and lower correlation coefficient in 
the Ghirlanda relation as compared to the Amati relation: 
TK,A = 0.67 ± 0.07 (5.6a) & tk,g = 0.66 ± 0.07 (5.5a) with 
aA = 0.18 & ac = 0.23. 

The same arguments in a yet stronger form hold for 
the comparison of the colhmation-corrected peak luminos- 



ity (Iz-y) with Ep^int correlation and the Lisa — Ep_i„t rela- 
tion (Ghirlanda et al. 2005b, hereafter G05b). No signifi- 
cant difference in the scatter of the two L-y — Ep^int & 
Lisa — Ep^int relations is observed when we consider 
the same sample of 16 GRBs with firm jet opening 
angle — reported in G04a — for both relations (0.15 
dex &: 0.16 dex respectively). 

How important are the <0.6% of certain outliers to the 
Ghirlanda relation? Unlike the case for the Amati relation, 
the method given by NP05a to determine the number of 
certain outliers, cannot be regarded as a rigorous test of the 
Ghirlanda relation, since it requires a knowledge of the jet 
opening angle distribution function as well as the redshift 
distribution of the bursts. Several attempts have been made 
so far to determine the distribution of jet opening angles 
(e.g. Frail et al. 2001; Norris 2001; Ghirlanda et al. 2005a; 
Guetta et al. 2005; Friedman & Bloom 2005), even given 
the sparse number of bursts with measured jet opening an- 
gles. However, as indicated by Band & Preece (2005), also by 
Perna, Sari & Frail (2003), a major problem with the current 
observed distribution of 9jst is that it is affected by another 
type of selection effect on its head and tail (i.e. very high and 
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very low 9 jet), which is related to the current limited abil- 
ity of observing very early and late breaks in the afterglow 
evolution of the bursts. Nevertheless, all obtained distribu- 
tion functions imply a range of Ojet < 40° for the jet open- 
ing angle with the peak of the distributions being around 
5° — 10°. Therefore, the extreme assumption that we made 
at the beginning of this section (i.e. ft = 1, corresponding to 
Ojet ~ 90°) in order to determine the number of outliers to 
the Ghirlanda relation appears to be unrealistically gener- 
ous. By this, we have assumed that the outflows of all bursts 
are essentially isotropic. Even if the relativistic outflows are 
not highly coUimated, some beaming is expected in most 
cases, since the energy channels mainly along the rotation 
axis of the inrushing material into the newly created black 
hole (Woosley, 1993). 

We therefore conclude that the use oi ft — 1 in the 
method given by NP05a, severely underestimates the num- 
ber of outliers to the Ghirlanda relation. In order to show 
how the apparent consistency of the bursts with Ghirlanda 
relation exacerbates using < 1 , we also find the Ghirlanda 
relation limit at > 3(j using the latest update of the relation 
(G07) and assuming Oj^t = 25° as the upper value of the 
jet opening angle - twice as large as the largest reported 
angle in the GRB sample of G07 - for which we find 7% in- 
consistency given the BATSE sample of LGRBs considered 
in this work. Another consistency check can be performed 
by estimating the beaming fractions (ft) of bright BATSE 
LGRBs via the relation proposed by Norris (2001) between 
the spectral lag of LGRBs and their beaming fractions. Even 
though we make generous assumptions in the calculation of 
the number of certain outliers, we find at least 34% of 1310 
bright BATSE LGRBs to be certain > 3cr outliers to the 
Ghirlanda relation (Eqn. 3 of G07). 

In addition, the redshift that maximizes the redshift- 
dependent term A{z, a, /?, ^, a) of NP05a's method in equa- 
tion Q is at Zmax > 10 for the Ghirlanda relation. This red- 
shift is much larger than the maximum detectable redshift 
by BATSE that Cohen & Piran (1995) report (2 = 2.llJ 7) 
assuming no evolution in the luminosity function of the 
bursts. Moreover, Norris & Gehrels (2008), have recently 
estimated the redshift distribution of SWIFT GRBs with 
unknown redshifts to be the same as the rest (1/3) of the 
SWIFT sample with measured redshifts resulting in an av- 
erage z ~ 2.1 for the whole sample of SWIFT bursts. Know- 
ing that BATSE LADs were on average 5 times less sensitive 
than BAT (Fenimore et al. 2004), we can use this as an up- 
per limit for the average redshift of BATSE GRBs which 
results in an A{z) that is a factor of two smaller than what 
was used to obtain the limits given in previous paragraph, 
leading to an increase in the number of certain outliers to 
the Ghirlanda relation. 

The consistency checks for this relation are, however, 
uncertain so long as the accurate unbiased redshift and jet 
angle distributions of the whole sample of BATSE bursts are 
not known. 



4 DISCUSSION 

The significant frequency, at over 19%, of > 3a outliers to 
the Amati relation found here is in contrast to the com- 
plete lack of outliers reported in G08 and the ~ 6% outlier 
frequency reported in N08. The presently reported outlier 
frequency is comparable to the ~ 25% outlier frequency re- 
ported by NP05a but still less than the ~ 88% outlier fre- 
quency reported by Band & Preece (2005). We find a similar 
outlier frequency to the Liso — Ep^i„t relation , > 8% at > 3cr, 
which is again in stark contrast to the findings of N08 with 
only ~ 0.2% at > 3<j outliers to this relation. All of these dis- 
crepancies can be attributable to two sources: either an old 
version of the Amati relation was being used, or a different 
subsample of BATSE LGRBs was being used. A third pos- 
sibility raised initially by Ghirlanda et al. (2005a), that the 
apparent high frequency of outliers found by the method 
given in NP05a, could be due to the assumption of a low 
scatter in the Amati relation, is rejected. Were it true, the 
practical uses of the Amati relation, as well as other cor- 
relations intimately connected with, most importantly the 
Ghirlanda relation, would be limited. 

Even if the Amati relation as given by G08 is exact, 
it is not a significantly more accurate estimator of redshifts 
than random, since the positions of the bursts relative to 
these relations appear to be indifferent to a wide range of 
z. Figure |9] (Left) shows the Eiso — Ep^i„t of a sample of 
LGRBs used by G08 to define the Amati relation together 
with their trajectories on this plane for a wide range of red- 
shifts (0.2 < 2: < 20). Inspection of the plot indicates that 
the scatter and correlation strength of the relation depends 
very weakly on the redshifts of the bursts. The same sample 
of bursts (black dots in the Left plot of Figure [9]l also define 
a strong correlation in the observer plane (green dots in the 
Right plot of Figure |9|. The trajectories are however, differ- 
ent at lower redshifts and deviate from the Amati relation in 
Eiso — Ep^int plane for all bursts. This could well explain why 
almost all outliers to Amati relation have redshifts z < 0.2. 

The Amati relation is created by apparently low disper- 
sion, highly correlated bivariate distribution of GRBs in the 
Sboi — Ep^obs plane. The low dispersion effect is itself cre- 
ated by the detection and selection effects on the faint and 
low energy edges. To show this more clearly, we derive the 
linear fits in the observer and rest frames to the sample of 
G08 bursts by excluding outliers in both planes as labeled 
in Figure [9] (except XRF 050416 A which is not an outlier in 
either of the planes and was labeled for another reason to 
be discussed below), for which we find, 

Log {Ep,obs) = 4.75 4- 0.55Lo<? (Sboi) , (10) 

Log {Ep,,„t) = -25.11 + 0.52Log {E,so) ■ (11) 

The variance between Stoi and Ep,int in both the ob- 
server and rest frames, as shown in the plots of Figure [9] 
has about the same scatter - 0.23 dex and 0.21 dex - 
with a slightly higher correlation coefficient being found in 
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Figure 9. Left: Plot of peak energy Ep^int versus the isotropically emitted energy Ei^o, with the G08 sample shown as black points. 
The apparent correlation of these points is the basis for the Amati relation. The dotted lines represent the trajectories of the LGRBs for 
different values of redshifts while holding all other GRB parameters fixed. The trajectories show the indifference of GRBs to redshift, 
indicating that the Amati relation is relatively insensitive to redshift. For any randomly given redshifts with z > 0.2, all G08 
LGRBs (black circles) reside within the 3cr consistency region of the Amati relation — the two blue dashed lines. In other words, there is 
no need to know the measured redshifts of LGRBs in G08 sample to define the tight Amati relation in the rest frame. 
One can almost always obtain a tight relation similar to the Amati relation in the rest frame by attributing each burst 
a random redshift. The head-down (empty red) triangle in the plot is XRF 0504f6A that is claimed by Sakamoto et al. (2006) to be 
a further confirmation of the Amati relation in the gap of this relation around the 30-80 KeV range. The solid green line is the redshift 
trajectory of this burst holding all other parameters fixed. The consistency of this event with the rest frame Amati relation, as noted in 
Figure [T] is not due to a physical origin, but merely a result of lying along the same region in the observer frame where LGRBs define a 
strong observer frame correlation [Right plot) due to selection effects. For any value of redshift z > 0.2, XRF 050416A is always consistent 
with the Amati relation at < 2a. Interestingly, all Amati relation outliers except GRB 050315 & 050826 have z < 0.2. The two empty 
red squares both represent GRB 060218 for two different reported spectral parameters by Campana et al. (2006) (shown as C06 in the 
figures) & B07. The same situation also holds for GRB 060614, as shown by the two red circles. The head-down red-filled triangle is GRB 
030329, consistent with the rest frame Amati relation, but an outlier to its observer frame counterpart. Right: Plot of Ep ^bs versus Stoi 
for the G08 sample of LGRBs used to define the Amati relation in the rest frame plane. The green colored data points represent the 
same LGRBs (black circles) of the Right plot. The labeled bursts are LGRBs that are outliers to either the observer frame or rest frame 
correlations (defined by the green and black circles in the corresponding plots), except XRF 050416A which is consistent with Amati 
relation in both planes. The consistency of XRF 050416A with the rest-frame Amati relation is due to residing in a region of observer 
plane that makes it also consistent with the observer-frame Amati relation. It is shown in the Left plot that the position of this burst in 
Eiso — Ep^int plane with respect to Amati relation is insensitive to its redshift and it will be consistent with the relation for any randomly 
given redshift z > 0.09. The head-up red triangle next to GRB 050315 is GRB 980425, a LGRB very well consistent with the observer 
plane Amati relation, but a far outlier to the rest frame relation. Both the observer and the rest frame Amati relations have 
comparable scatters (0.23 dex &: 0.21 dex respectively) and correlation coefficients that are not significantly different 
from each other (rK.obs = 0.59 ± 0.05, 7.5(t & TK,rest = 0.65 ± 0.04, 8a respectively). Except XRF 050416A, GRB 050826 & GRB 
050315, all other labeled GRBs have a redshift z < 0.2. Any burst with z > 0.2 that is consistent with the observer frame 
relation will always be consistent with the rest frame relation, while a burst with z < 0.2, consistent with either of the 
two relations will always be an outlier to the other. 



the rest frame of the bursts: TK,obs ~ 0.59 ± 0.05, 7.5a & 
TK,rest = 0.65 ± 0.04, 8(T respectively. This slight improve- 
ment, however, is statistically marginal, undermining a po- 
tential physical origin to the Amati relation. According to 
F-test, there is only la {p — 0.285) weak evidence of a sig- 
nificant difference between the variances of the observer and 
rest frame Amati relations. 

In order to show how redshifting of the parameters on 
both sides of the relation Eqn. ( 10 1 can boost the existing 



ulation by giving each LGRB used to construct the relation 



( 10 1 a random redshift taken from the sample over a large 



correlation in the observer plane, we run a Monte Carlo sim- 



number of iterations taking into account of the limited en- 
ergy budget of LGRBs, assuming = 0.27, f^A = 0.73 & 
Ho ~ 72 Km/sMpc. Averaging over all iterations we find 
that the Kendall's rank correlation coefficient of the ob- 
served sample is generally enhanced by ~ 0.05, that is about 
the same as the correlation coefficients difference (~ 0.06) 
obtained for the sample of bursts with real measured red- 
shifts. Also, 53% of the iterations result in a median scatter 
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in the rest frame distribution that is smaller than the median 
scatter of the observer frame. 

Inspection of Figures [t] &[9](ie/i plot), however, indi- 
cates a lack of bursts on the lower right sides of the observer 
& rest frame planes of Sboi — Ep,obs & Eiso — Ep^int ■ In other 
words, GRBs that are both bright and soft appear to be 
rare. This could have a physical origin and would likely be 
unaffected by detection threshold limits. In our estimation, 
the existence of a probable physically-based lower 
limit for the hardness (Ep^obs) of GRBs as a function 
of fluence, together with the selection effects on the 
hard dim side, have possibly led to the creation of 
the low-dispersion Amati relation, and similar Ep^int- 
based correlations such as the Ghirlanda relation, in 



the rest frame of the bursts (Figure 10 1 



It is interesting to wonder ^why almost all of the out- 
liers to the Amati relation appear to be sub- energetic, lying 
at redshifts of z < 0.2?' Inspection of Figure |9] {Left) indi- 
cates a likely reason: Ep^int is insensitive to redshift when 
the redshift is small. Specifically, any burst that lies within 
3(7 limits of this relation in the observer plane with a redshift 
z > 0.2, will always be consistent with the Amati relation 
(Eqn. 11 1 in the rest frame at < 3a regardless of the any 



quent outliers at lower isotropic energies due to a 
general insensitivity of Ep^i„t to redshift at jz < 0.2. 



5 CONCLUDING REMARKS 

Throughout the presented analyses, we investigated the ef- 
fects of BATSE Large Area Detectors' threshold limita- 
tions on triggering GRBs over a wide range of spectral 
peak energies (Ep^obs)- In order to compare the distribu- 
tion of BATSE GRBs with these triggering thresholds on 
the plane of bolometric fluence vs. peak energy, we relied on 
the hardness — Ep^obs correlation reported by Shahmoradi & 
Nemiroff (2010a) to map the hardness vs. bolometric fluence 
(Sboi - HRh) plane of 1900 BATSE GRBs - for which accu- 
rate continuous light curve data as well as fluence in all en- 
ergy channels were available - into the plane of Sboi —Ep^obs- 
The results of the simulation indicate that the distribu- 
tion of GRBs on Sboi — Ep^obs plane is affected by the BATSE 
trigger thresholds on the regions where dim hard bursts re- 
side. This is simply due to the fact that BATSE speciflcally 
and GRB detectors generally are photon counters rather 
than bolometers. Therefore, for a given bolometric peak flux 



possible redshift that the burst might have. On the other (e.g. in units of ergs cm ^ s ^), the harder bursts would 



hand, any burst consistent with relation ( 10 1 at < 3(t with 
a redshift z < 0.2 would almost always be a certain outlier 
to Amati relation at > Sa in the rest frame, or in reverse, 
any burst with z < 0.2 that is consistent with the Amati 
relation, would almost always be a certain outlier to rela- 



tion ( 10 1 at > 3(7. Therefore, the apparent deviation from 



the Amati relation observed for sub-luminous bursts is un- 
physical and created by the random redshifts of the bursts. 
It is also of no surprise that the newly detected bursts are 
generally consistent with this relation, so long as they are 
detected and selected in the region of consistency with the 
observer frame Amati relation with a z > 0.2. 

As an example, we consider XRF 050416A that is re- 
ported by Sakamoto et al. (2006) to be a further "conflrma- 
tion" of the inclusiveness of the Amati relation in the gap 
around the 30-80 KeV range. The redshift trajectory of this 
burst is shown in Figure [o] (ie/t plot), holding all other pa- 
rameters fixed (the green solid line in the graph). Given any 
random redshifts {z > 0.1), this burst is always consistent 
with Amati relation at < 3a, indicating a nonphysical origin 
for its consistency. 

It is interesting that there are no nearby {z < 0.2) bright 
outliers to Amati relation which is possibly due to the lumi- 
nosity distribution and the redshift evolution of luminosity 
function of the bursts (e.g. Kocevsky & Liang 2006; Lloyd- 
Ronning et al. 2002) as well as possible yet unknown effects 
of jet opening angles of the bursts on their detections. In 
other words, the lack of bright bursts detected in the 
very nearby universe results in the apparent tight- 
ness of the Amati relation at high isotropic energies 
{Eiso ^ 10^^ [ergs]) as well as a large scatter with fre- 



have less photon count rates than the softer GRBs and con- 
sequently, less chances of detection. Although the signifi- 
cance of the BATSE triggering thresholds on the bivariate 
distribution of LGRBs in Sboi—Ep^obs plane has yet to be de- 
termined (Shahmoradi & Nemiroff 20106), the current sam- 
ple of LGRBs detected by BATSE provides evidence against 
the proposed low-dispersion Amati & Liso — Ep^int relations 
(Yonetoku et al. 2004; Amati 2006): 

About 20% of the BATSE sample of LGRBs considered 
here appear inconsistent with the Amati relation at > 3a 
based on the traditional (Tgo > 3 [sec]) or fuzzy clustering 



classification (Shahmoradi 2010) of BATSE LGRBs (jj2.2l. 
The fraction of outliers is likely more than the ratio obtained 
for reasons of sample inclusion and statistical methodology. 
The current sample would likely include more outliers were 
some GRBs themselves not excluded by detection thresholds 
of BATSE LADs. In addition, the method of NP05a used 
herein can only set a lower limit to the number of the out- 
liers. This is bolstered, knowing that none of the currently 
known outliers to the Amati relation with firmly measured 
redshifts, except GRB 050826, could be flagged as certain 
outliers at > 3a by this method. 

Similar conclusions also hold for the simulations of 
BATSE trigger thresholds on the plane of Pboi — Ep^oba {{ 2.3 1. 
Using the same method as applied for Sboi — Ep^ots plot of 
BATSE GRBs, we can set a lower limit on the inconsistency 
of BATSE LGRBs with Liso - Ep,,nt relation. We find that 
at least 8% of the bright BATSE LGBRs are certain outliers 
to this relation as given recently by G09 at > 3a significance 
level. 

It is notable that the results of our simulations are in 
contrast with the findings of previous authors, in particular 
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Figure 10. A plot of 1-second peak flux P1024 versus Ep ^bs for 
1900 BATSE GRBs. The main plot contains only LGRBs clas- 
sified via a fuzzy clustering algorithm described in Shahmoradi 
(2010), while the inset plot contains only SGRBs. The dashed 
line represents the nominal BATSE trigger threshold on 1024ms 
timescale. The lack of very soft bright bursts is visible in the top 
left of both plots which seems to have a physical origin, while 
there are apparently no limits for the positions of the bursts on 
the dim hard side of the plots except BATSE trigger threshold. 



N08 & G08, where they find that BATSE trigger thresh- 
old, as well as spectral analysis limits on the two planes of 
Sboi — Ep^obs & Pboi — Ep^obs have possibly little or no ef- 
fects on the distributions of BATSE LGRBs on these two 
planes. The reason for the discrepancies should be sought 
in the values of the limiting parameters that they use in 
their simulations, such the average ratio of fiuence to peak 
flux of the bursts (FPR) and the nominal durations of the 
bursts considered therein, all taken from a small fraction of 
detected, spectrally analyzed BATSE LGRBs. In this sense, 



their analyses suffer from a circular logic problem {\2.2\. In 
addition, the use spectral models with fixed photon indices 
in their simulations, results in a severe underestimation of 



the selection effects {{ 2.3 1. The strong evolution of the peak- 
energy (Ep^obs) in the light curves of the bursts (e.g. K06, 
Ryde 1999; Crider et al. 1999; Band 1997; Crider et al. 1997; 
Liang & Kargatis 1996; Ford et al. 1995) is another factor 
that is overlooked in the simulations of G08 & NOB. We have 
worked to make the simulations presented here free from the 
above mentioned deficiencies. 

Another strong argument that favors an unphysical ori- 
gin for the Amati, Ghirlanda and possibly other 3-parameter 
relations, such as the empirical Liang-Zhang relation, comes 
from inter-comparisons between the proposed relations. Pre- 
vious authors have reported a significant scatter reduction 
in transforming the Amati relation to these 3-parameter re- 
lations, specifically the Ghirlanda relation. However, consid- 
ering the same sample for both relations that are being com- 



pared to each other, we find that the scatter reduction and 
correlation improvements are insignificant (!|3|. Therefore, 
in order to have a meaningful comparison of any two rela- 
tions with each other, in particular the Amati & Ghirlanda 
relations, it is important to consider the same data set for 
both relations. 

It is also noteworthy that the sample of LGRBs con- 
sidered by G08 to construct the tight Amati relation also 
shows a strong correlation in the observer frame, with a 
scatter comparable to the dispersion in the rest frame Am- 
ati relation, differing by only 0.02 dex. This indicates that 
the tightness of the Amati relation is only a ghost of the 
tight correlation of this sample of LGRBs in the observer 
frame, reinforced by redshifting the spectral parameters of 
the bursts from the observer to the rest frame plane. For 
any random redshifts that these bursts might have, the rest 
frame Amati relation is on average always tighter than the 
Amati relation in the observer frame (jQ. Moreover, the ap- 
parent frequent inconsistencies of the sub-luminous LGRBs 
with the Amati relation, appear to have no physical origin 
and can be attributed purely to redshifting the spectral pa- 
rameters of the bursts that mainly reside on a narrow strip in 
the observer frame, by a redshift z < 0.2 (*j4]& Figure [7] |9| . 

In sum, the Amati relation as proposed by Amati 
(2002), Amati (2006) & Ghirlanda et al. (2008) appears to 
be greatly affected by complex selection effects in triggering, 
spectral analyses & redshift measurements of LGRBs on the 
dim side of the Sboi — Ep^obs plane. The lack of LGRBs on 
the soft bright side of the Stoi — Ep^obs plane might possibly 
retain an underlying physical origin. Nevertheless, the 
practical use of Ep^obs as a standard candle is questioned, 
as its detector convolutions likely compromise its use as a 
discerning probe for cosmological models. 



This work could not have been accomplished without 
the vast time and efforts spent by many workers over the 
past decade, in particular BATSE team, including the 
designers, builders, and analysts for the burst detectors on 
board the Gamma Ray Observatory who have accumulated 
and analyzed the observations and summarized them in 
BATSE GRB catalogs. In particular, we acknowledge sev- 
eral useful communications with David Band, and dedicate 
this paper to his memory. 
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